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Abstract 

This paper deals with ratio-dependent predator-prey systems with delay. We will inves- 
tigate under what conditions delay cannot cause instability in higher dimension. We give an 
example when delay causes instability. 



Key words and phrases: predator-prey system, functional response, sign stability, ratio dependence, 
delay 
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1 Introduction 

Let us consider the following ratio-dependent ecological system, in which n different predator 
species (the i-th predator quantities at time t are denoted by yi(t), i = 1, 2, . . . , n respectively) are 
competing for a single prey species (the quantity of prey at time t is denoted by x(t)): 

x = rxg(x 

i=l 



i=i 

Vi = ViVi (^) - diVi, i = l,2,...,n 



x . 

where dot means differentiation with respect to time t. We assume that the per capita growth 
rate of prey in absence of predators is rg(x, K) where r is a positive constant (in fact the maximal 
growth rate of prey), K > is the carrying capacity of environment with respect to the prey, 
the function g satisfies some natural conditions, see the details in [6]. For example one of these 
conditions is the following: 

(K-x)g(x,K)>0, x>0, x^K. (2) 

Such a function g is the so called logistic growth rate of prey 

g{x,K) = l-^. (3) 

We assume further that the death rate dj > of predator % is constant and the per capita birth 
rate of the same predator is pi (—), where the function pi also satisfies some natural conditions, 



x 

see also in [6]. 

In that paper we have already investigated the system with the Michaelis-Menten or Holling 
type functional response in case of ratio-dependence: 
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and with the ratio-dependent Ivlev functional response: 

Pi (^> a i) = m i (l - e - wj , (5) 

where parameter cij is the so called "half- saturation constant", namely in the case where Pi is a 
bounded function for fixed aj > 0, rrii = sup pi(x,yi,cii) is the "maximal birth rate" of the z-th 

x,yi>0 

predator. That means, if the functional response is a Holling-type without ratio-dependence then 
aj means the quantity of prey at which the birth rate of predator i is half of its supremum. In case 
of a ratio-dependent Holling model etj means a proportion of prey to predator at which the birth 
rate is half of its supremum. In case of an Ivlev model the meaning of a\ is similar to the earlier, 
see the details in [6]. (To save space we did not write out the dependence on aj in (1).) For the 
survival of predator % it is, clearly, necessary that the maximal birth rate be larger, than the death 
rate: 

m> di- (6) 

This will be assumed in the sequel. Finally, we assume that the presence of predators decreases 
the growth rate of prey by the amount equal to the birth rate of the respective predator. 



2 Model with delay 

We get a more realistic model if we take into account that the predators' growth rate at present 
depend on past quantities of prey and therefore a continuous weight (or density) function / is 
introduced whose role is to weight moments of the past. Function / satisfies the requirements: 

POO 

f(s)>0, s€(0,oo); / f(s)ds = l, (7) 

Jo 
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and x(t) is replaced in the growth rate of predator i by its weighted average over the past: 

q(t):= [ x(r)f(t-r)dr. 



This means that the time average of prey quantity over the past has the same fading influence on 
the present growth rates of different predators. The simplest choice is f(s) = ae~ as , with a > 0. 
This function satisfies the condition (7) and now 

q(t) = [ x{T)ae- a{t - T) dT. (9) 



We call this choice of / exponentially fading memory, see in [2], [7]; later in [4]. (Since / is the 
probability density of an exponentially distributed random variable, the probabilistic interpretation 
is obvious.) The smaller a > is the longer is the time interval in the past in which the values of 
x are taken into account, i.e. - is the "measure of the influence of the past". It is easy to see that 
with this special delay, system (1) is equivalent to the following system of ordinary differential 
equations: 

n 

x = rxg(x,K) - ^yiPi (— 



i=l 



Vi = ViPi ( — ) - diUi, i = 1, 2, . . . , n 



(10) 

1 

q = a(x — q) 

where function Pi(-) can be (4), (5) or any kind of general ratio-dependent functional response if 
we replace x(t) by the time average q(t) of prey quantity over the past. Similar systems have been 
studied by many authors in the two-dimensional case, specially in [1], and also with diffusion in 
[8]. In [1] the functional response was of the simplest Holling-type one without ratio-dependence 
and in [8] the functional response was of the Michaelis-Menten-type with ratio-dependence and 
also with diffusion. Our aim in this paper is to study the effect of exponentially fading memory in 
case of a general ratio-dependent functional response with more than one different predators. 
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The qualitative behaviour of (1) was studied in [6], where it has been supposed that there exists 
an equilibrium point E*(x*,yl, . . . ,y*) in the positive orthant, where x*, and y* are the solutions 
of the following equations: 

n 

rxg(x,K) = s ^d i y u pi(^j=d u i = l,...,n. (11) 

8=1 

Note that x* > if and only if K > x* because of (2). 
The coefficient matrix of the system (1) linearized at E* is: 



A 



an 



?,* n'* ( Vn ~ 1 ' 



where 



an 



Pi 



-dx-y\v'th 



-dx-ylv'i^ 




o 










o 



Un-lPn-l }* 



rg{x 



K)+rx*g' x (x*,K)-J2y;P?(-^ 



2>i 



i=l 



Pi - 

x 



VnPn x * 



(12) 

(13) 
(14) 



An n x n matrix A 



is said to be stable if each of its eigenvalues has a negative real part. 



The following definition can be found in [5]: 

Definition 2.1. An n x n matrix A = [ay] is called sign-stable if each matrix A of the 
sign-pattern as A (signay = signa^ for all is stable. 



same 
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It was proven in [6] the following: 



Theorem 2.1. // 



and 



a u < 0, 



Pi =Pi [ — ) < 0, i = l,...,n, 



-di-y*p'*— = -di-ytp'i ( ^- ) — < 0, i = 1, ...,n 

X \ X / X 



(15) 
(16) 

(17) 



then matrix (12) is sign-stable, thus, E* is an asymptotically stable equilibrium point of system 
(!)■ 

Now, let us suppose that there exists a positive equilibrium point E*(x*, y{, . . . , y*) of system 
(1), then with the definition q* := x* and E^(x*, y{, . . . , y*, g*) we get an equilibrium point of (10) 
in the positive orthant. And again x* > if and only if K > x*. 
The coefficient matrix of system (10) linearized at E* d is: 



an 





di-ytPx^ -d 2 



ylp';^ 



ylp'i^ 



* /* 1 

y 2 P2^ 



~d n VnPn x * 





y*v'* — 







i/Spf(-£) 



y*nPn( 



a: 











-a 



, (Vi 

Pi ( - 

V x 



d V i (f) 



4t) 



where an is given by (13) and again p'* = p\ y~ 
We note that (18) can not be sign-stable because its graph have cycles. (See in [5].) 
Let us restrict the number of predators to two. 
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2.1 One prey two predators with delay 

Let us consider system (10) in case of n — 2. We suppose that (15), (16), (17) hold for i — 1,2. 
In this special case the entries of matrix A d are a n < 0, a 22 ,a 33 < 0, a 12 ,ai 3 < 0, a 24 ,a 34 > 0, 
a 4i = a > 0, a 44 = —a < 0. This means that A d has the following sign pattern: 



-/0 - - 

- + 

0- + 

a —a 



The characteristic polynomial of a matrix with the same sign pattern as (19) is: 



(19) 



.D(A) = A 4 + a 3 A 3 + a 2 A 2 + a x A + a 



(20) 



with 

«3 = -an - ct22 - ct33 + a, 
a 2 = ana 22 + ana 33 + a 22 a 33 - a(a n + a 22 + a 33 ), 
ai = -ana 22 a 33 + a(ana 22 + ana 33 + a 22 a 33 ) - a(ai 2 a 24 + ai 3 a 34 ), 
a = det A d = a(-ana 22 a 33 + a 22 ai 3 a 34 + a 33 ai 2 a 24 ). 
It is known that the necessary condition of stability of the polynomial -D(A) is > 0, % = 0, 1, 2, 3. 

Lemma 2.1. // yLj has the same sign pattern as (19) then the above necessary conditions of 
stability are satisfied for all a > 0. 



Proof. It is an elementary calculation to prove a, > 0, i = 0, 1, 2, 3, for all a > 0. 



□ 
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Sufficient condition of stability of matrix in this case is: 

03(0102 — O0O3) — a\ > (21) 

See for example Theorem 1.4.8 in [3]. It leads to a very complicated formula. In order to check 
this we used Wolfram Mathematica 6.0. http:/ /www.wolfram.com. We got: 
H(a) = 03(0102 — O0O3) — a\ — 

( 2 2 2 2 o 2 2 

— O22O11 — 033011 — a 22 Oii — OggOn + Oi20240n — 2022033011 + 013034011 — 0220 33 + O12O22O24 — 2 2 a 33 + 

O13O33O34) a 3 

+ (022ofi + a33afi + 2a2 2 Oii + 20330ii —012024011+4022033011—013034011 +02 2 0n +033011 +4a22o| 3 on — 

2 3 2 2 o 2 2 

O12O22O24O11 + 4o 22 0330n + 012024033011 + 013022034011 — 013033034011 + 0220 33 — a 12 a 24 ~r 2o 22 3 3 + 

222 2 3 2 2o 

012024033 — Oi 3 034 — O12O22O24 + O22O33 + O12O22O24O33 + 013022034 — 013033034 — 2012013024034 + 

\ 2 

013022033034; a 

+ (— a^Oii —o^Oii — 2022033011 — a| 2 Oii — o^o^ — 4022033011 +O12O22O24O11 — 4a| 2 0330ii +013033034011 — 

o3/i22 2 2 o 3 2 

2O22O33O11 — 4a 2 2033an — O12O24O33O11 + O12O22O24O11 — 2a 2 20330n — O12O22O24O33O11 — 013022034011 + 

2 23 332 23 2\ 

a 13 a 33 a 340ll —013022033034011 —022033 — 012024033—022033 — 012022024033—013022034 — 013022033034; a 

+Oll02 2 o| 3 + O11O22O33 + Oiia2 2 o| 3 + 20ii02 2 o| 3 + Oiia 2 2o| 3 + Oii02 2 033 + Oii02 2 033 

Lemma 2.2. // matrix (18) in case of n = 2 has a pure imaginary eigenvalue then in (21) the 
expression at left hand side is equal to zero. 

Proof. If we substitute ju, j 2 = —1, u ^ into (20) we get u 2 = ^ and 03(0102—0003)— Oi = 0. □ 

As we can see by result of Wolfram Mathematica 6.0 the left hand side of condition (21) has 
the following form depending on a: 

H(a) = A 3 a 3 + A 2 a 2 + A x a + A (22) 
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Lemma 2.3. If A^ has the same sign pattern as (19) and an < then A 3 , A > 0. 

Proof. The proof is complete by elementary calculations. □ 

Lemma 2.3 means that the function H(a) given by (22) is positive, and monotone increasing 
or decreasing depending on A\ > or A\ < 0, respectively; and has a convex or concave down 
shape if A 2 > or A 2 < 0, respectively; at a = 0. 



la lb lc 




Figure 1: The value of A\ and of A 2 is positive 



Figures 1, 2, 3, 4 show that there are several cases when delay does not destabilize the system 
for any a, for example if A 2 > 0, A\ > 0, and the cases when H(a) has a single real root only. 
Furthermore, if a increases through a limit, namely if - is small, "measure of the influence of the 
past" is small then the system (10) has a locally asymptotically stable equilibrium point E* d . This 
situation corresponds to our expectation and it is similar as it was in the 2-dimensional case, see 
in [1]. 

Now we can formulate our main result. We will give appropriate conditions that can easily be 
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2a 2b 2c 




Figure 2: The value of A\ is positive and of A 2 is negative 



3a 3b 




Figure 3: The value of A\ is negative and of A 2 is positive 
checked in order to satisfy A 2 > 0, A\ > 0. 

Theorem 2.2. If matrix A^ given by (18) in case of n = 2 satisfies conditions (15), (16), (17) for 
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4a 4b 




Figure 4: The value of A\ is negative and of A2 is negative 

i — 1, 2 (has the same sign pattern as (19)) and the following two conditions also hold 

a n > a 33 > -^13034, (23) 

°11 > a 22 > -^12^24 (24) 

then Ad is stable and is an asymptotically stable equilibrium point of the delayed system (10) 
in case of n = 2 for any a > 0. 

Proof. Under the conditions of the theorem we can decompose the expression of A\ into the 
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following positive terms: 

A\ = (aj 2 + Ctl2'3'24)( — «33 — O44O33 — QII.O22O33) + (O33 + 043034) ( — «22 — a ll a 22 — O44O22O33) 

+ ( a ii _ a3 3 )( a 22ai2024) + (an _ 032X033043034) 

+ / 3 2 23i 2 r)3 a 2 2 3 

^ — a n a 22 — a, ll a 22 ' a ll a 22 a 12 a '24 — ^ a n a 22 a 33 — ^ a n a 22 a '33 — a ll°22 a '33 

32/12 2q22 23 3 i 2 \ 

— 011033 — ^0^-^0,22033 — /CI11G22O33 — O11O33 — O44O22O33 ~T~ O44O33O43O34 ) 

> 

and similarly for the expression of : 

A2 = (d\ 2 + O12O24X011O33 + O22O33 + Ogg — O12O24) + (O33 + 013034) (O44O22 + O22O33 + O22 — 013034) 

+ ( — 011012024 — 0^043034 — 2012024043034) 

/ 3 2 2 3 3 2 

+ (044022 + 2044022 + 0440 2 2 — 044022042024 + O44O33 + 4a 11 022033 + 

2 2 2 2 3 \ 

30440 2 2033 + 2044033 + 3044022033 + 044033 — 044033043034) 

> (a 22 + O42O24X044O33 + 022033 + 033 — O42O24) + (033 + 043034) (044022 + O22O33 + O22 — 043O34) 

+ ( — 033042024 — 0^043034 — 204 2 a 2 4043a34) 

/ 3 2 2 3 3 2 

+ (0440 22 + 20440 22 + 0440 22 — 044a 22 04 2 a 2 4 + O44O33 + 40440 22 033 + 

2 2 2 2 3 \ 

3a44a 22 a33 + 2044033 + 3a44a 22 a 3 3 + 044033 — 044033043034) 

= ( a 22 + a 4 2 2 4) (044033 + 22 033 + 033 — O42O24) + (033 + 043034) (0440 22 + 22 033 + a 22 — 043034) 

+ (~o 12 a 2 4(a33 + 043034) — 043034(032 + 042024)) 

/ 3 2 2 3 3 2 

+ (O44O22 + 2044022 + 0440 2 2 — 044022042024 + 044033 + 4044022033 + 

2 2 2 2 3 \ 

30440 2 2033 + 2044033 + 3044 022 33 + 044033 — 044033043034) 

> 

□ 
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This theorem means that in case of a sign-stable interaction matrix (12) there are many cases 
when delay does not destabilize the system. By Theorem 2.1, if an < (given by (13)) and if 
conditions (16), (17) are also satisfied then (12) is sign-stable. This is the two-dimensional situation 
modeled by Farkas and Cavani in [1] when the equilibrium point lies on the descending branch of 
the prey nullcline. That is the case when E* lies outside the Allee-effect zone - here the effect 
of overcrowding is already felt. Any further increase in prey quantity must be counterbalanced 
by a decrease in predator quantity, see in [4]. On the other hand, in the Allee-effect zone prey is 
scarce and an increase in prey quantity is beneficial for the growth rate of prey, see in [4]. Let us 
introduce the vector 



Vector (25) has two rows Fi and F 2 in the two-dimensional case. Suppose that any predator 
quantity growth will decrease the growth rate of prey, namely F[ y < 0. Some typical reasonable 
forms of Fi(x,yi) = zero isoclines applicable to most species in case of ratio-dependence are 
shown in Figure 5. We can see that F[ > 0, thus an > in the Allee-effect zone modelled by the 
increasing branch of the function in the third graph. 

In case of our model we keep this meaning of the Allee-effect zone, and we say we are outside 
of the Allee-effect zone if in order to keep the prey growth rate zero the increase of prey can be 
counterbalanced by the decrease of the whole quantities of the different predators. Let us consider 
the higher dimensional cases. Now the function F given by (25) has n+1 rows Fi, i — 1, 2, . . . , n+1. 
Suppose that any predator quantity growth will decrease the growth rate of prey, namely F[ < 0, 




i=l 





(25) 
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0.02 0.04 0.06 0.08 0.10 q.00 0.02 0.04 0.06 0.08 0.10 02 04 006 00s 010 



Figure 5: Typical nullclines of prey in case of 2-dimension 



% = 1,2, ... ,n. In the three dimensional case a typical onion-like prey zero isocline surface of 
Fi(x,yi,y 2 ) = is shown in Figure 2.4.2 in [4] page 44 without ratio-dependence. Inside the 
onion-like surface F\ > while outside F\ < 0. Function F is increasing as we cross the surface 
inwards and therefore its gradient points inward. Therefore if the equilibrium point is on the 
eastern hemisphere of this onion then F[ < 0, thus, an < and on the western hemisphere of the 
onion F[ > 0, thus, an > and we can see that F[ > 0, thus an > in the Allee-effect zone. 
The onion is similar to this in case of ratio-dependence shown in Figures 6, 7, 8. 

If F[ < (namely yi is predator of x) then a n > holds also in higher dimension in the Allee- 

n 

effect zone. To see this, let us consider F\{x,yi, . . . ,y n ) = rxg{x,K) — ^^yiPi ^— J and surface 

i=i x 

F\{x, yi, . . . , y n ) = 0, which is the prey zero isocline surface. Let be E 1 = (x 1 , y\, . . . , yl), E 2 = 
(x 2 , y\, . . . , y 2 n ) two different points in the Allee-effect zone on the prey isocline surface, where 
x 1 < x 2 , y\ <y 2 , ,i = l,...,n. 
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0.010 




0.10 

Figure 6: Typical zero-cline of prey in case of r = 3 in 3-dimensions (r = 3, K = 0.1, m\ = 
16, a\ = 4, m 2 = 18, a 2 = 2) 

= F 1 (x 2 ,y 2 1 ,...,y 2 n )-F 1 (x 1 ,yl,...,y 1 n ) 
= {F^x 2 , y x 2 , . . . , yl) - F^x 2 , yl y 2 2 ..., y 2 n )} 

r Z7i f 2 1 2 2\ T7/21122 2m 

+ l^iO ,2/i,2/ 2 --->2/J -Fi{x ,y 1 ,y 2 ,y 3 ,y 4 ,...,y n )\ 

+ {^i0 2 , 2/i\ 2/2> 2/3> 2/4 > • • • > 2/n) - Fi(x 2 , yl, yl, y\,yl,...,y 2 n )} + ... 

+ {Fi{x 2 , yl, yl, yl,..., ?/*_i! yl) - F i( x2 , vl, yl,---, vl-i, vl)} 

+ {F^x 2 , yl, yl, yl,..., y l n ) - F,{x l , yl, . . . , £)}. 
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0.10 

Figure 7: Typical zero-cline of prey in case of r = 7 in 3-dimensions (r = 7, K — 0.1, mi = 

16, ai = 4, m 2 = 18, Gt2 = 2) 

Expressions in the brackets are negative except the last bracket because of F' y . < 0, thus F x > 
must hold. 

It is reasonable to say that E* lies outside the Allee-effect zone if an < and E* lies in the 
Allee-effect zone if an > 0. 

Remark 2.1. Theorem 2.2 means, if E* lies outside the Allee-effect zone then delay does not 
change the stability behaviour of the system in this special case. 

This remark is a direct generalization of Case 1 of [1] on page 226. 
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0.10 

Figure 8: Typical zero-cline of prey in case of r = 10 in 3-dimensions (r = 10, K — 0.1, m\ = 
16, a\ = 4, m 2 = 18, a2 = 2) 

The meaning of conditions (23), (24) are the following: 

Conditions a\ x > a^ 3 , > a| 2 mean that intraspecific competition in prey species is greater than 
intraspecific competition in predators species. 

The meaning of conditions > —013034, a| 2 > —0-12024 is in connection with the phenomenon of 
their consume strategy, namely do they try to ensure their survival by having a relatively high or 
low growth rate and are able or not to raise their offspring on a scarce supply of food. We will 
discuss this very interesting meaning of conditions (23), (24) in case of (3) and (4) or (5) in the 
following section. 
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2.2 Strategies 

The condition an < can be ensured by a relative high intrinsic growth rate r of prey. This means 
that there is enough food for predators in order to reproduce well. If this fact is valid in a long 
term then we expect even more that a predator species has an advantage that need more food and 
has a high growth rate. The parameter aj > is the half saturation constant of predator i. This 
means that when the quantity of prey reaches value a« then the per capita birth rate of predator 
% reaches half of the maximal birth rate, as one can see in case of a simple Holling model where 
Pi(x, di) = Wj^q^, rrii is "the maximal birth rate" of the i-th predator, and Pi(di, aj) = In case 
of ratio-dependent models parameter a« has a similar meaning, namely the greater ctj is the more 
food is needed for predator i. To see this let us consider the ratio-dependent Holling function, 
given by (4). In this case at a fixed value of yi, Pi(x, yi, ai) = ^ if x = Oiy». Similarly in case of the 
ratio- dependent Ivlev function, given by (5) at a fixed value of yi, Pi{x, yi, a*) = ^ if x — a^yi In 2. 
Thus, a predator with a big half saturation constant can be considered as an r-strategist and with 
a lower one as a K-strategist. See in [6], [4]. Thus, we expect that the parameters cannot be 
arbitrary small, because the mentioned effect is stronger in that case when the time average of prey 
quantity over the past has the same influence on the present growth rates of different predators. 
The following theorems express this situation. 

Theorem 2.3. Let matrix be given by (18) in case ofn = 2 satisfying conditions (15), (16), (17) 
for i = 1,2 (i.e. has the same sign pattern as (19)) and the function g, Pi are given by (3), 
(4), respectively. If Oi > 1 fori = 1,2 then conditions (23), (24) are satisfied. 

Proof. Calculate a| 3 > — ai3a 34 , a\ 2 > — ai 2 a.24 by substituting (3), (4) and the statement follows. 

□ 
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Theorem 2.4. Let matrix be given by (18) in case ofn = 2 satisfying conditions (15), (16), (17) 
for i — 1,2 (i.e. Ad has the same sign pattern as (19)) and the function g, pi are given by (3), 
(5), respectively. If ai > | fori = 1,2 then conditions (23), (24) are satisfied. 

Proof. Calculate a| 3 > —013034, a\ 2 > — 012O24 by substituting (3), (5) we get: 

rrii mi mt—di 



{> (ln^V) 2 • [b) 



Let us denote x = ™\ , x > 1. Thus, 

mi— dj' ' 

1- i - Mill 

where lim^i+o ai(x) = ~ and ai(x) is monotone decreasing for x > 1 because its derivative is: 

/ \/ ((In x) 2 -2^+2+2 In a;) . . . . 

aj(a;) = (i^ip anc ^ the numerator is negative because it is zero if x — 1 and the 

derivative of ((lnx) 2 — 2x + 2 + 2 lnx) is negative for x < 1. Thus, the maximum of the righthand 
side of (26) is equal to | and theorem holds. □ 

The meaning of Theorems 2.3, 2.4 corresponds to our expectation, namely in case of delayed 
models the advantage of the r-strategist can be seen over the K-strategist. This advantage is 
greater in case of a ratio-dependent Holling model than in case of a ratio-dependent Ivlev model. 
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2.3 One prey, n predators with delay 



Now let the number of predators n be an arbitrary positive integer and let us consider system (1} 
with its coefficient matrix given by (12). Let us denote the entries of (12) by a^-, thus 



A 



an 

021 

a 3 i 



a a2 ... 

a 2 2 
a 33 



Oln 






On-1,1 



a n-l,n-l 



(27) 



a ni ... a nn 
If we modify system (1) with delay we get system (10) which, after linearization has the coefficient 
matrix given by (18). We have seen that (18) can be obtained from the entries of A as follows: 



A n 



an O12 Ol3 

a 22 
a 33 



Oln 

a 2 i 

03i 



nn Q"nl 

—a 



(28) 



... a 
a ... 

Theorem 2.5. Let matrix Ad be given by (18) for arbitrary positive integer n, and suppose it 
satisfies conditions (16) and (17) for all i = 1,2, ... ,n; and let an < 0. If a is small enough or 
large enough then A^ is stable, and E* d is an asymptotically stable equilibrium state of the delayed 
system (10). 

Proof. Let us consider the characteristic polynomial T>(\) := det(y4^ — XE) of (28). Let us denote 



column i of matrix Ad by Cj 



1,2, ... ,n) and let us make the following column operations: 
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first c x =>- Ci + c n+ i, then c„ +1 =>- c n+1 — ci. Now we get 

an — A fli2 

Q>21 a 22 ~ A ... 

«3i a 33 - A 

det(A d - A£) 



din -(an - A) 





Ojnn. A 



-a: 



Let us make the following partition of this determinant: 



det(A d - XE) = det 



A-XE 



-(an - A) 




-A 



... 



-a 



det 



(29) 



C D 



(30) 



Applying the Schur theorem [9, Theorem 3.1.1] we get: 



det(A d - XE) = det (A - XE) det(A d - XE\A - XE), 



where (A d — XE\A— XE) is the Schur-complement of A— XE in A d — XE, namely (A d — XE\A — XE) 
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D — C(A — XE) 1 B and suppose that A is not an eigenvalue of A. 
(A d - XE\A - XE) = D-C{A-XE)- l B 



-a-[ -X 



]{A- XE) 



-i 



-(an - A) 




—a — A(an — X)A U , 



where : = 



det(A-AE) ( fl22 _ A ) ' ' ' 

det(A d - AE|A - XE) — -a — X 



a nn - A), thus, 

(an - A) (a nn - A) 



det(A-XE) 



We get the following relation (true for all A e C) 



det(A d - XE) = -a det(A - XE) - X(a n - A) (a nn - X) 

n 

= (-l)(adet(A- XE)+XY[(a u -X)). 



i=l 



(31) 



(32) 



Now we prove that the coefficients of this polynomial have the same sign, using the fact that A 
being sign stable, hence the coefficients of det(A — XE) have the same sign. Let us denote the 
coefficients of det(A — XE) by a i: namely: 



det(A - XE) = (-X) n + a n _i(-A) n_1 + • • • + a . 
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Thus, 



det{A d - XE) 



(-l){a(-A) n + aa n _i(-A) n_1 H h aa 



+A((-A) n + (a n + • • • + a nn )(-X) 



n-1 



+ (ana 22 H h a„_i ri -ia rin )(-A) 



H h (a n a 2 2 




(-A) n+1 + (a n + ■ ■ • + a 



■nn 



a)(-A) 



71 



+(ana 2 2 + • • ■ + a 

n— In— 1 



a 



aa n _i)(-A) 



n-1 



'mi 



H h (ana 2 2 



a 



"iirt 



aai)(— A) — aao- 



Since det(A — XE) is a stable polynomial, hence if n is even, then a 2 fc is positive, and a 2 fc + i is 
negative for all k. Thus, the coefficients with even indices of det(Ad — XE) are negative, and those 
with odd indices are positive, and all the coefficients of (A)- 7 (j = 0, 1, . . . , n + 1) in det(A d — XE) 
are negative. 

For the case of n odd we can repeat the above proof. Thus the necessary condition of stability 
of the polynomial det(Ad — XE) holds. 

This means that if det(Ad — XE) is not a stable polynomial then it has to have a pair of complex 
conjugate roots with nonnegative real part. 

Now let us consider the case when a is very large. Then the eigenvalues of det(Ad — XE) 
are close to the eigenvalues of A and there is a remaining root with an unknown sign. But this 
root should also be a negative real number, because it has no pair to be a member of a complex 
conjugate pair, and because the coefficients of the characteristic polynomial are positive. Thus, 
for sufficiently large a ^> the matrix A d is stable. 

If a is very small then the eigenvalues of det(Ad — XE) are close to the roots of A IlT=i( aii — ^) = 
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0. It has n negative real roots and one more root left with an unknown sign. And again, this should 
be a negative real number, because it has no pair to be a member of a complex conjugate pair, and 
because the coefficients of the characteristic polynomial are positive. Thus, for sufficiently small 



The meaning of this theorem is the following. If a is small then the measure of the influence 
of the past is large. In this case the equilibrium point E^ is locally asymptotically stable. 

If a is large then the measure of the influence of the past is small, the system's behaviour is 
close to the behaviour of the system without delay, of which the equilibrium E* was stable. Thus, 
the results correspond to our expectations. But all these are true outside the Allee-effect zone, 
where the stability is stronger than inside. 

2.4 Numerical examples 

Example 2.1. Let us consider a three dimensional Holling type ratio-dependent model with 
delay, namely g is given by (3) and Pi is given by (4). Let the constants be given as follows: 
m\ = 16, rri2 = 18, d\ = 8, c?2 = 12, a\ — 4, a.2 = 2, K = 0.1. The equilibrium point of 
the system depending on r is E* = (0.1(1 — |), ^(1 — -), ^(1 — -)) . In this case the interaction 
matrix of the system without delay is given by: 



q^O the matrix is stable. This completes the proof of the theorem. 



□ 



8-r -4 -8 



.4 



1 



-4 



(33) 



1 







4 



The characteristic polynomial of A is: D(\) = (—4 — A) (A 2 + (r — 4) A + 4(r — 5)) . This is a stable 
polynomial for r > 5 and A is sign stable for r > 8 . 
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The equilibrium point of the delayed system depending on r is 

E d * = ( 0.1(1 - -), 1(1 - -), 1(1 - -), 0.1(1 - -)) . 
\ r 40 r 40 r r J 

The coefficient matrix of the delayed system linearized at Ed* is 



— r 


-4 


-8 








-4 





1 








-4 


1 


a 








—a 



(34) 



The characteristic polynomial of Ad is: D d (\) = (—4 — A) ((8 — r — A) (—4 — A) (—a — A) — 12a) . 

Let us check conditions (23), (24). It is easy to see that in case of r > 12 these are satisfied. 
The conditions of Theorem 2.2 hold, Ed* is asymptotically stable. Time evolution of the species 
are shown on the left side of Fig. 9, whereas the right side shows the corresponding trajectory 
together with the equilibrium point. The form of (22) with r = 13 is shown in Fig. 10. This 
corresponds to Fig. 1, case lc. It is easy to see that the equilibrium point of the delay system 
remains asymptotically stable for any a > 0. We note that in this case the equilibrium point is 
outside the Allee-effect zone, see Fig. 8. 

If 12 > r > 5 then conditions (23), (24) are not valid, and there are such cases when Ed* is 
stable and there are cases when it is unstable. Time evolution of the species are shown on the 
left side of Fig. 11, whereas the right side shows the corresponding trajectory together with the 
equilibrium point. The form of (22) with r = 7 is shown in Fig. 12. It is easy to see that there are 
values of a for which H(a) < 0, thus, the equilibrium point of the delay system is unstable, and 
also values for which H(a) > 0, thus, the equilibrium point of the delay system is asymptotically 
stable. We note that in this case the equilibrium point is inside the Allee-effect zone, see Fig. 7. 
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Figure 9: Left: Time evolution of the species in case of r = 13, a = 1. Right: The trajectory 
tends to the asymptotically stable equilibrium point, (x is red, q is green, yi is dashed blue, ?/ 2 is 
yellow.) 




Figure 10: The function (22) with r = 13 



Ratio-dcpcndcnt predator- prey systems with memory 



28 




Figure 11: Left: Time evolution of the species in case of r = 7, a — 1. Right: The trajectory 
leaves the neighbourhood of the unstable equilibrium point, (x is red, q is green, y\ is dashed blue, 
y 2 is yellow.) 




Figure 12: The function (22) with r = 7 



Of course this study is not complete. There are many interesting trajectories, periodic orbits, 
see e. g. Fig. 13, 14. 
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Figure 13: Left: Time evolution of the species in case of r = 7, a — 1. The solution seems to be 
periodic at first sight. (The reason of this phenomenon may also be numerical errors.) Right: The 
corresponding trajectory, {x is red, q is green, y 1 is dashed blue, y 2 is yellow.) 




Figure 14: Left: Seemingly time periodic evolution of the species in case of r = 8, a = 0.2. Right: 
The corresponding periodic orbit, {x is red, q is green, yi is dashed blue, y 2 is yellow.) 
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The interested reader can experiment with the parameters and initial conditions of the model 
using the Mathematica program on the page http:/ /www.math.bme.hu/~jtothhtml#kktj. E. g. it 
is also interesting how the trajectories change if we reduce aj. In case of r < 5 there is no positive 
equilibrium point E%. 

The mentioned program produces figures like Fig. 15. In case of an Ivlev model similar 




Figure 15: Snapshot of manipulation, r = 8 and a = 0.25. 
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situations may occur. 
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